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Abstract 

We study ion condensation onto a patterned surface of alternating charges. The competition be- 
tween self-energy and ion-surface interactions leads to the formation of ionic crystalline structures 
at low temperatures. We consider different arrangements of underlying ionic crystals, including 
single ion adsorption, as well as the formation of dipoles at the interface between charged do- 
mains. Molecular dynamic simulation illustrates existence of single and mixed phases. Our results 
contribute to understanding pattern recognition, and molecular separation and synthesis near pat- 
terned surfaces. 
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The interaction of molecules with surfaces or interfaces is a fundamental problem in 
modern science. Advances in molecular design and maturity of nanoscale microscopy paved 
the way in a number of novel technological achievements, such as the ability;^ to_ produce 
nano-structured surfaces , and complex self-assembled aggregates 
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riodic and spatially organized surfaces can be used as templates for fabricating nanostruc- 



tures, magnetic storage media, nanowires 



or substrates for control of cell adhesion and 



growth 



10|. From a scientific point of view, charged surfaces and their interactions with 



ionic environments are particularly interesting due to the long-range nature of electrostatic 
interactions and its relative independence from molecular structure. The break down of the 
mean field Poison-Boltzmann approximation to describe adsorption of ions to charged sur- 



faces 



ll| motivated many studies. Simple models such as a single component plasma near a 



neutralizing background revealed two-dimensional crystallization of adsorbed ions onto ho- 
mogeneously charged surfaces That is, the electrostatic repulsion between condensed 
ions results in lateral correlations. 

In this letter we study ion condensation on a planar surface with periodic domains of 
alternating charge. In particular, we consider striped patterns of width A and constant 
surface charge number density a (Fig. [11(a)). The half-space above the plane is filled by a 
solution of positive and negative ions with valency Z, radius a and bulk number density 
p° = p\ = p°_. We consider a symmetric system with ions of equal valency. At low 
temperatures, the homogeneously charged regions of the surface attract oppositely charged 
ions. The long-range electrostatic interaction induces per iodic arrangement of condensed 
ions within the plane in analogy with Wigner crystals 13|]. The competition between ion- 
ion and ion-surface electrostatic interactions, together with the geometrical restrictions set 
by the surface patterns determine the structure of the ionic crystal. When the ion-ion 
interaction dominates, we find the formation of dipolar pairs along the interfaces between the 
surface patterns (Fig.[T](b)). On the other hand, when the ion-surface interaction dominates, 
the ions localize where the electrostatic potential is strongest, at the center of the stripes 
(Fig. [D^c)). It is the goal of this letter to formulate an approach that captures the essential 
physics of the phenomenon and gives a reasonable description of ion condensation near 
patterned surfaces. 

We model the adsorption and the interaction of the ions between themselves and with the 
surface as follows. The free energy of the absorbed ions is purely electrostatic and includes 
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contributions from interactions between the ions and surface as well as between ions. The 
free energy of the ions in the bulk is dominated by their translational entropy. We note that 
the bulk ions are screened from the surface by the condensed ions within the Gouy- Chapman 



length [l^, lo = l/27TaZ£. This allows us to consider the layer of adsorbed ions separately 
form the bulk. Equilibrium between the two regions is achieved when the chemical potentials 
of condensed ions, /ic, and bulk ions, fi, are equal fic = yW- Equivalently, this condition can 
be stated as the requirement that the excess energy of the condensed ions F, per unit area 
A, be a minimum with respect to the condensed ion density and arrangement. We write 



this effective free energy as: 



lF(p,/i) fiZ^M , 

-p{^ +v + e-fi], (1) 
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where i = /AneokBT is Bjerrum length, p is a surface number density of condensed ions, 
M is the Madelung constant of a crystal of condensed charged units, n is a number of ions 
per unit, v is the energy of ion-surface interaction, and e is the self-energy of the condensed 
units. The condensed unit can be single ion (n = 1) with zero self energy, e = 0, dipole 
composed of two opposite charges {n = 2), with e = —Z'^i/a, or more complex structure 
like a chain or a cluster of ions. 

The minimum of the free energy is determined by evaluating and comparing the values 
of expression Eq. [T] for different ion densities and lattice symmetries. The key step in these 
evaluations is the determination of the Madelung constant for the basic arrangements of 
ions (Fig. [It^a,b)). This constant is the effective potential experienced by a charged unit due 
to the presence of a lattice of similar units. For the system of dipoles (Fig. [11(b)), we have 



ij ij J 



ij ij 

where i and j are lattice indices, p is a dipole momentum of magnitude p = 1/a, Eij is 
the electric field at the position fij = {h{j + i/2), Xi). The sum over integers i and j runs 
from — oo to +oo and the prime over summation indicates the avoidance of singularity at 
the vector = (0,0). Due to the constant width of the stripe, the separation distance 
between condensed ions along the stripe can be written as h = n^/p+X, where is the 
number of positively charged ions in the Wigner crystal cell. For the case of ion absorption 
to the center of the stripe (Fig. [11(c)), we obtain: 



where ipij is the potential at position fij = {jh,\i). These sums can be calculated usin, 



Ewald summations [14], following, for example, the methods developed by Crandall 15] 



For the dipolar case, which does not pose convergence problems, faster techniques are avail- 



able 
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171]. Besides these two basic arrangements, we have also calculated Madelung 
constants for lattices with a relative translation along the stripe for adjacent lines of ions. 
We find that both lattices shown in Fig. [T]^b,c) have the lowest energies. 
The electrostatic potential near the striped surface (Fig{T](a)) is given by 

oo 

dx / dy^ iJi{x,y,Ay,z) (4) 

-oo -A/2 '=-°° 

where e is an electron charge (negative), a is a surface number charge density, A is the width 
of a stripe and 

(-1)* 

ilji(x,y, Ay, z) = — , (5) 

^ ' 'x^ + {y-Ay + i\f + z^ ^ ' 



where Ay e [—A/2, A/2] determines the shift from the center of the stripe along the y— axis. 



Using the identity 



171 ] the electrostatic potential can be calculated exactly: 
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^sin2(7rZ/2) /vr/fA - 2A'y)\ .i. 



P V 2A , 
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The magnitude of the potential decays exponentially and has sinusoidal profile with ex- 

tremum in the center of the stripes. Thus, the energy of the single charge in the center of 

the stripe (Fig. [U^c)) is Vc = — 47rZe^craA/3 and the energy of the dipole at the interface 

(Fig. [1Kb)) is Vd = -2nZe^aa^. 

In order to investigate the phase behavior of the system first we determine the conditions 

for ion adsorption into each of the two states. The adsorption of ions into a regular structure 

starts when the bulk chemical potential is lower than the critical value implicitly determined 

from the equation: 

dp \A ksT ) ^„ 

when we consider the free energy as a function of the chemical potential. It is equivalent 
to appearance of the first non-trivial extremum of the free energy. This conditions (Eq|7]) 
can be solved exactly for two types of condensed ionic lattices. In the limit of small surface 



charge density, the adsorption starts from the formation of the dipoles along the boundary 
of the surface patterns (Fig. [T](b)), due to the weak interaction with the surface. The critical 
chemical potential for the system of condensed dipoles reads 

It is important to note that it is A independent. With further increase in the surface charge 
density, the energy of the electrostatic attraction to the surface overcomes the energy of the 
dipole formation and the ionic crystal of condensed ions changes its structure. The new 
ionic crystal (Fig. [H^c)) has a simple rectangular lattice and is formed from the single ions 
located in the minima of the surface electrostatic potential. The critical chemical potential 
in that case is 

,1 . _« ( ^-^-^^f°^Pl ) (9) 

The transition from one structure to another occurs at critical surface charge density 

O'er ^ 3 / A - 2a log [2] \ 

Z 2Txa\ V 4A - 3a y ' ^ ' 

when the energy of both systems are equal, or equivalently, both have the same chemical 
potential. The critical surface charge density determines the beginning (cxcr, /icr) of the 
coexistence line yUt^, which separates regimes for condensed ions forming ionic crystals of 
dipoles (Fig. [I](b)) from the one of the ionic crystals of single charge (Fig. [H^c)). 

Once one of the condensation conditions /i < is satisfied, the associated ions form 
structures with well defined densities p. In these conditions, the density p of the condensed 
ions can be calculated from 



for the corresponding lattice symmetry. The preferred symmetry of the ionic crystal is 
determined by direct comparison of the two energy minima achieved by each of the structures 
considered. We note that at the transition between structures, when both energies are equal, 
the associated density is different, and the transition between phases is first order. 

The stability of phases A and B is determined as follows. Phase A is composed of dipolar 
units, and in general, the addition of one more dipole does not reduce the energy. However, 
the introduction of an isolated ion at the center of an oppositely charged stripe might 
indeed reduce the energy, signaling the fact that phase A is no longer the minimum-energy 



conformation. Similarly, we consider the addition to a background of stripe centered ions, 
phase B, of a single dipole at the interface between two oppositely charged stripes. When 
these perturbations reduce the energy of the phase with least energy (between A and B), 
the actual minimum of that state corresponds to a more complex structure. It is likely that 
the actual minimum structure is a combination of both striped-centered ions and dipoles, 
denoted here as a mixed phase. For example, for the dipole phase A, the instability criteria 
reads: 

SF iZ'^ , , 

Mm + t^c - /i < 0, (12) 



where Vc is the energy of the ion-surface interaction, and 

oo 

= iTxapA + levraAV^ ^ lK^l}?p^ cos (vr/) (13) 

z=i 

is the energy of a single ion at the center of the stripe in the presence of the dipolar array 
(Fig. [It^b)), and is a modified Bessel function of first order. The spinodal boundary 
(when the system is first unstable), appears at values of the chemical potential //^ for which 
bF = 0. The spinodal boundary for the stripe-centered case is labelled /x^. 

We illustrate these results for the cases A/a = 10 and Zi = 1. Figure [2](a) shows the 
dependence of the critical chemical potential vs. surface charge density cr for systems of both 
symmetries (Eqs. [H]and[9]). The condition of equal chemical potential between phases A and 
B determines the coexistence boundary line (fitr)- Figure [2]^a) also shows the location of the 
spinodal curves for the dipole /x^, and the centered /i^ states. Along the coexistence line 
between phases A and B, the system undergoes a first order phase transition associated with 
the jump in the density of condensed ions. Figure [2](b) shows the dependence of the density 
for both systems along the coexistence line. Figure [3](a) summarizes results in a schematic 
phase diagram in terms of the surface charge density a/Z and the chemical potential fi/Z"^ 
normalized by the valency Z of the ions. 

Figure [2](a) shows two basic features of the boundaries of stability for the dipolar and 
centered structures. Note that the spinodal line for the stripe centered state, system B, 
lies approximately at a constant chemical potential. At lower values of the potential, it is 
simply more convenient to add single ions at the center of the stripe. Reciprocally, the limit 
of stability of the dipolar state, system A, appears at near constant value of the surface 
charge. This surface charge density determines the strength of the electric field at the 
interface between stripes, and sets the energy gain for a deposited dipole. For strong enough 



values of the charge density, the energy of the stripe centered state can be lowered by adding 
dipoles at the boundaries. Large values of both, chemical potential and charge density, lead 
therefore to the complex mixed state in the upper right corner of our phase diagram. 

We explore finite temperature effects in the phase diagram with molecular dynamics 
simulations. Condensed ions form dipoles along the interfaces and single ions lie in the 
center (Figl3](b,c)). The strong correlation between ions, is preserved even at relatively high 
temperature T = CSfc^T; thermal fluctuations do not destroy the long range order. In our 
simulations we have found, in addition to dipole and single charge phases, chain clusters 
and paired dipoles suggesting rather complex phase behavior previously found mostly in the 



bulk 
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Our exploration of the low temperature limit of ion condensation onto a structured surface 
shows that there are many different controllable behaviors that a fixed surface pattern can 
generate. It seems plausible to use the properties of condensed ions to create complex, 
tunable structures at interfaces that exhibit a variety of potential applications. For example, 
electrostatic charge on membrane surface can be used to attract free vesicles or membrane- 
coated charged microbeads 20| to employ intermembrane adhesion and lipids exchange 
with substrate-supported patterned membrane. Another example, a new more effective and 
selective method, compared to previous DNA fractionation techniques, was demonstrated 



2l|. 



using streaming dielectrophoresis of DNA on the Si surface with micro-scaled Au strips 
Control of the abundance of the charged species can lead to dramatic transitions, that 
can be used as catalytic steps of complex chemical reactions, as steps in the formation of 
secondary self-assembled structures superimposed on a simpler pattern. Though our analysis 
is focused only on two basic behaviors, stripe-centered, and dipolar arrangements, we show 
the existence of more complicated, mixed, phases. We expect, in fact, a cascade of transitions 
between all these complex phases. 

This work is supported by NSF grant numbers DMR-0414446 and DMR-0076097. The 
authors thank G. Vernizzi for helpful and stimulating discussions. 
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FIG. 1: (a) The schematic representation of the system. The basic Wigner lattices formed by 
condensed (b) ions forming dipoles and (c) single ions. 
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FIG. 2: (a) The dependence of the eritical chemieal potential {^cr)-, eoexistence {jitr) and spinodal 
{^Jisp) hnes vs. surface charge density a/Z. (b) The dependence of the density vs. the surface charge 
density for both Wigner lattices formed by dipoles {pa) and single ions [pb) along the coexistence 
line. 
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FIG. 3: (a) The schematic phase diagram imphed by results of calculation, where lines (1) and 
(2) correspond to and /i^, line (3) is a coexistence line ^tr-, line (4) and (5) are spinodal lines 
^fp and /x^ correspondingly, (b, c) Snapshots from molecular simulation of mixed and single ion 
phases at T = O.bksT. 
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